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Abstract. Recent observations by XMM-Newton detected rotational pulsations in the total brightness and spectrum of several 
neutron stars. To properly interpret the data, accurate modeling of neutron star emission is necessary. Detailed analysis of the 
shape and strength of the rotational variations allows a measurement of the surface composition and magnetic field, as well 
as constrains the nuclear equation of state. We discuss our models of the spectra and light curves of two of the most observed 
neutron stars, RX J1856.5— 3754 and IE 1207.4—5209, and discuss some implications of our results and the direction of 
future work. 
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INTRODUCTION 

Thermal radiation from the surface of neutron stars (NSs) 
can provide invaluable information on the physical prop- 
erties and evolution of NSs. NS properties, such as the 
mass M and radius R, in turn depend on the poorly 
constrained physics of the stellar interior, such as the 
nuclear equation of state (EOS) and quark and super- 
fluid/superconducting properties at supra-nuclear densi- 
ties. Many NSs are also known to possess strong mag- 
netic fields {B ~ 10'2 - IQi^ G), with some well above 
the quantum critical value (B » Bq = 4.4 x lO''' G). 

The observed thermal radiation originates in a thin at- 
mospheric layer (with scale height ^ 1 cm) that cov- 
ers the stellar surface. To properly interpret the obser- 
vations of NS surface emission and to provide accurate 
constraints on their physical properties, it is important 
to understand in detail the radiative behavior of NS at- 
mospheres in the presence of strong magnetic fields (see 
fl, 2, 3, 4, 5], for more detailed references on observa- 
tions and on previous works in NS atmosphere model- 
ing). The properties of the atmosphere, such as the chem- 
ical composition, EOS, and radiative opacities, directly 
determine the characteristics of the observed spectrum. 
While the surface composition of the NS is unknown, 
a great simplification arises due to the efficient gravita- 
tional separation of light and heavy elements |6]. A pure 
hydrogen atmosphere is expected even if a small amount 
of accretion/fallback occurs after NS formation; the total 
mass of hydrogen needed to form an optically thick at- 
mosphere can be less than ^ lO'^ g. On the other hand, a 
heavy element atmosphere may be possible if no accre- 
tion takes place. 



The strong magnetic fields present in NS atmospheres 
significantly increase the binding energies of atoms, 
molecules, and other bound states (see |7], for a review). 
Abundances of these bound states can be appreciable in 
the atmospheres of cold NSs (i.e., those with surface tem- 
perature T < 10^ K; |8, 9]). In addition, the presence 
of a magnetic field causes emission to be anisotropic 
and polarized; this must be taken into account when de- 
veloping radiative transfer codes. The most comprehen- 
sive early studies of magnetic NS atmospheres focused 
on a fully ionized hydrogen plasma and moderate field 
strengths (B - lO'^ - lO'^ G; iH [U [H [11). These 
models are expected to be valid only for relatively high 
temperatures. More recently, atmosphere models in the 
ultra-strong field (B > lO'"* G) and relevant ternperature 
regimes have been presented ( ifTil [Tsi 0, [IB. [H; see 
also [18, 19], for early work), and all of these rely on 
the assumption of a fully ionized hydrogen composition. 
Magnetized non-hydrogen atmospheres have been stud- 
ied by 1 10, 20], but because of the complexity of the 
atomic physics, the models were necessarily crude (see 
i2ll], for more details). Only recently has self-consistent 
atmosphere models ll22il23ll2lll using the latest EOS and 
opacities for partially ionized hydrogen Il24l25[l and mid- 
Z elements I 26ll27ll been constructed. 

The atmosphere models discussed above only describe 
emission from a local patch of the stellar surface. By 
taking into account surface magnetic field B and tem- 
perature T distributions, we can construct more phys- 
ically correct models of emission from NSs. However, 
these spectra from the whole NS surface are necessarily 
model-dependent, as the B and T distributions are un- 
known. Nevertheless, detailed comparisons of the mod- 



els with rotation phase-resolved observations is a pow- 
erful tool to study NSs, e.g., spectral features that vary 
with phase are essential to disentangling magnetic field 
effects from other parameters and to probe the magnetic 
field geometry on the surface of the star Indeed there 
have been recent works attempting to fit magnetic atmo- 
sphere spectra to observations of NSs (see |28], and ref- 
erences therein). Here we describe some of the details 
and observational applications of our work. 

MODEL FOR NEUTRON STAR 
SURFACE EMISSION 

In order to construct models of emission from the entire 
NS surface, we first build atmosphere models for a given 
effective temperature Teff and magnetic field strength B 
and direction 0^ relative to the surface normal; these 
(local) models describe a patch of the NS surface. In 
the presence of magnetic fields B > lO'^ G, radiation 
propagates in two photon polarization modes (see, e.g., 
II29II ): therefore, the atmosphere models are obtained by 
solving the radiative transfer equations for the two cou- 
pled polarization modes (see |2, 22, 30, 23], for details 
on the construction of the atmosphere models). In addi- 
tion to B and Tgff, the atmosphere models have a depen- 
dence, through hydrostatic balance, on the surface grav- 
ity g [= (1 +Zg)GM/R^], where the gravitational redshift 
Zg is given by (1 +Zf,) = (1 - 2GM / Rc^y^l^; however, 
the resulting spectra do not vary significantly using dif- 
ferent values of g around 2 x lO'* cm s^^ 1 1]. 

Next, the entire NS surface is divided into regions with 
different B and Teff. Relatively simple surface distribu- 
tions of B and Teff are adopted: we assume the surface is 
symmetric (in B and T^s) about the magnetic equator and 
divide the hemisphere into several magnetic colatitudinal 
regions. An example parametrization is given in Table[T] 
note that the magnetic field distribution is roughly dipo- 
lar. Emission from any point within a colatitudinal region 
is given by the atmosphere model for that region. 



TABLE 1. Neutron Star Surface Parametrization 
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Finally, the spectra from the entire NS surface is com- 
puted by summing over the emission from the differ- 
ent regions for a given rotation phase (see llsil l28ll . for 
details). The observed emission depends on two angles 
(a, C,): a is the angle between the rotation and magnetic 
axes and C, is the angle between the rotation axis and the 



direction to the observer. We also account for the bend- 
ing of the path of light due to gravity, which causes more 
of the NS surface to be visible (135° and 115° as com- 
pared to 90° without Hght-bending for M = 1 AMq and 
R= \Q km and 14 km, respectively; see ll32il33ll ). 
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FIGURE 1. Light curves for different geometries (a, ^): 
class I with (20°,30°), class II with (30° ,60°), class III with 
(60°, 80°), and class IV with (10°,80°). a is the angle between 
the spin and magnetic axes, and C, is the angle between the 
spin axis and the line-of-sight. The four classes are defined in 
1 33]. The solid lines are for the magnetic model described in 
the text [dashed line is for (50°, 50°)], while the dotted lines 
are analytic light curves (scaled arbitrarily in amplitude) for 
isotropic emission from two antipodal hot spots (see | 33]). 

Figure[T]shows the light curves, for various geometries 
(a, 0, of the NS model using the parametrization given 
in Table [T] Zg = 0.2, and g = 1.1 x lO'^ cm s^l We 
also plot the analytic light curves from for isotropic 
emission from two antipodal hot spots (see ll34[[35[] . for 
examples of pulse profiles from non-magnetic hydrogen 
atmosphere hot spots). The classification scheme (for 
isotropically-emitting hot caps) is defined in ll33ll : (I) 
only the primary cap is visible, and the pulse profile is 
purely sinusoidal with a single peak, (II) the opposite cap 
is seen around pulse minimum due to relativistic light- 
bending, which reduces the strength of the modulation, 
(III) the primary cap is not seen during a segment of the 
rotation, and (IV) both spots are seen at all phases and 
thus there is no modulation. 

Several important features are evident from a com- 
parison of magnetic atmosphere emission to that of 
isotropic emission. The angular-dependence of the radi- 
ation (or beam pattern) manifests as a narrow "pencil- 
beam" along the direction of the magnetic field and a 
broad "fan-beam" at intermediate angles (see lfT2i [T^, 
for beam patterns and spectra at various ©5). As dis- 
cussed in il2ll . the pencil-beam is the result of the lower 
opacity at angles < {E /EbY/^, where Eg — heB/nieC = 
11.6(B/10'2 G) keV is the electron cyclotron energy; 
the width of the pencil-beam is thus {E/EbY^^, and 



the radiation is more strongly beamed at higher mag- 
netic fields. This narrow beam is seen in the (50°,50°)- 
light curve plotted in FigurefT] which is the only instance 
shown that has the observer's line-of-sight exactly cross- 
ing the magnetic cap and coinciding with the peak of the 
isotropic emission. Also evident is the fan-beam (most 
obvious in the light curves of classes II and III), which 
occur on either side of the magnetic cap and can increase 
the number of light-curve peaks. Finally, the anisotropic 
beam pattern (combined with the surface temperature 
variation) can produce an apparent phase shift compared 
to isotropic emission and modulation when an isotropic 
beam pattern suggests none (c.f. class IV). 

RX J1856.5-3754 

RX J1856.5 — 3754 is one of the brightest, nearby isolated 
NSs ll36l] . and considerable observational resources have 
been devoted to its study. Recently, XMM-Newton ob- 
servations uncovered pulsations from RX J1856.5— 3754 
with a period of 7 s and an upper limit on the period 
derivative IHIMl. The EPIC-pn and MOS Ught curves 
are shown in Figure |2l Several characteristics of the ob- 
served light curves suggest possible values of a and 
(1) the 1.6%±0.2% amplitude of the pulsations, (2) a sin- 
gle peak (or two peaks close in phase, as may be the case 
for the pn light curve) per rotation, and (3) no significant 
energy-dependence for the single observation and small 
energy-dependence for all XMM-Newton data ll28lf37ll . 
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FIGURE 2. Energy-integrated (0.15-1.2 keV) light curves 
of RX J1856.5-3754. Histograms are the XMM-Newton EPIC- 
pn (top) and EPIC-MOS (bottom) observations (see Fig. 1 of 
137,11 ). Solid and dotted lines are the models with (3°,40°) and 
(7°,60°), respectively, along with the fit deviations [i.e., (data- 
model)/a] in the corresponding lower panels. Two rotation 
periods are shown for clarity. 



In previous work [30], we found that our models 
of a magnetic, partially ionized hydrogen atmosphere 
matches well the entire spectrum, from X-rays to opti- 
cal, of RX J1856.5— 3754 [best-fit parameters: gravita- 
tional redshift Zg ~ 0.22, magnetic field B « 4 x lO'^ G, 
and radius R" w 17 km, where R°° = R{ \ + Zg)\ see Fig- 
ure[3l. With the discovery of rotational modulation of the 
X-ray emission 1 37], we use the light curves predicted by 
our model (described in the previous sectionVto constrain 
the geometry (a, of RX J1856.5-3754 Q. We find 
angles of < 6° and « 20 — 45° (example fits ares shown 
in Figure |2|i. This indicates that either the rotation and 
magnetic axes are closely aligned or we are essentially 
seeing down the spin axis of the NS (see also I39ll ). 




FIGURE 3. Spectrum of RX J1856.5-3754 from optical 
to X-ray wavelengths. Data points are observations taken from 
various sources; error bars are IcT uncertainties. The solid line 
is the absorbed (and redshifted by Zg = 0.22) atmosphere model 
spectrum with S = 4 x lO'^ G, T^ff = 5.3 x lO'' K, and R°° = 
17 km; the dashed line is the unabsorbed atmosphere model 
spectrum. The dash-dotted line is the (absorbed) blackbody fit 
to the X-ray spectrum with = 5 km. Note that our atmo- 
sphere model underpredicts the optical flux by 15%; however, 
observational and model uncertainties here is ~ 20%. 

From the results of our modeling of the phase- 
resolved observations, we can better determine 
where RX J 1856. 5— 3754 belongs in the broader 
context of NS studies. For example. Figure |4] shows 
the constraints placed on the mass-radius relation- 
ship for NSs. Our results imply a relatively stiff 
but standard nuclear EOS (see, e.g., [40]). Also if 
RX J 1856. 5— 3754 is losing rotational energy by 
magnetic dipole radiation, the rate of spindown is 
given by dP/dt = 10"'^ s s-i(B/10'2 Gf{P/\ s)-i = 
5 x 10^'^ ss^'; this is well below the upper limit 
dP/dt < 1.9 X s s"^ obtained by jQ and illus- 

trated in Figure |5] Note that RX J1856.5-3754 has 
not been detected in the radio fisl HI, IH], though its 
location in P P space (see Figure |5]l is below the 
theoretical death line for radio pulsars. 



FIGURE 4. Constraints on NS mass M and radius R from 
fitting the observations of RX J1856.5— 3754. The upper left 
regions are excluded by general relativity and causality. The 
upper right and lower left regions exceed the ~ 30% uncertainty 
in R°°, which is dominated by the uncertainty in the distance 
(see 1 36]). The dotted lines indicate 3<J uncertainty in Zg', note 
that this is just the uncertainty from the fit and does not include 
systematic uncertainties in the data and model. 

IE 1207.4-5209 

IE 1207.4-5209 is a NS (with spin period = 0.424 s) 
in the center of supernova remnant G 296.5+10.0. Its X- 
ray spectrum is remarkable in that it shows two broad 
absorption features at 0.7 and 1.4 keV ll46i l47i ksl 
imiloi- Even more surprising is that these features show 
greater phase-variability than the continuum spectrum 
IdTI I52I1 . Proposed models involving ion cyclotron or 
atomic lines from a light element atmosphere at B ~ 
10^4 G 0111 seem unlikely due to weakening of line 
strengths by vacuum resonance effects |4, 54, 17] or low 
abundance of the ionization states possibly responsible 
for the observed lines ll27ll . Recent timing analysis also 
imply B < 3.3 x 10^' G (by assuming vacuum dipole 
braking); this suggests the spectral features are electron 
cyclotron lines |55]. Alternatively, lIsH, IztIi proposed 
a mid-Z element atmosphere at B ^ lO'^ G. We are 
studying this last case by building atmosphere models 
composed of mid-Z elements (see [21], for details) in an 
attempt to fit the phase-resolved observations. 

As an illustration of our models. Figure |6] shows 
the phase-resolved model oxygen atmosphere spectra, 
light curve, and pulse fraction [= (Cmax - Cmin) / (Cmax + 
Cmin), where C is the count spectrum]. We assume here 
that Zg = 0.2, (a, Q = (10°, 20°), and the NS surface is 
parametrized by r^ff = 4 x 10*^ K and [B{10^^G),&b]= 
(1.2,0°) for magnetic colatitudes 0-5° and (1,30°) for 
colatitudes 6-20°. Note that our models roughly match 
the continuum shape and line locations and strengths of 



-3 -2 -1 1 

log P (s) 

FIGURE 5. NS period derivative P as a function of 
spin period P. Lines indicate constant magnetic field [B = 
3.2 X 10''' (/>/>) G]. Dots are pulsars whose P and P are 
obtained from http://www.atnf.csiro.au/research/pulsar/psrcat/ 
(see iiUl l. We highlight a particular class of NSs, for which 
RX J1856.5— 3754 is a member: sources with a measured B 
from their X-ray spectrum (open squares) and the same sources 
with a measured P (solid squares) or upper limit to P (tri- 
angles) (see |42], and references therein). The star indicates 
RX J1856.5-3754 from our measured B. 



the phase-averaged spectrum of IE 1207.4—5209. More 
detailed analysis will be reported elsewhere. 
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FIGURE 6. Top: Model atmosphere spectra (convolved with 
the XMM-Newton EPIC-pn energy resolution) at different rota- 
tion phases. Inset: Energy-integrated light curve. Bottom: Pulse 
fraction as a function of energy. The dashed horizontal lines in- 
dicate the observed pulse fractions over the given energy range 
for IE 1207.4-5209 (see 1521). 

In summary, we discussed briefly our continuing work 
on understanding the emission process of NS surfaces. 



as well as detailed comparisons of our models to phase- 
resolved observations of the NSs RX J1856.5-3754 and 
IE 1207.4—5209. Analyses of other sources is ongoing. 
In the near future, we are integrating our model spectra 
into XSPEC for community use IIStIi . 
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